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ABSTRACT 

A general purpose micromechanics analysis that discretely models the yam architecture 
within the textile repeating unit cell, was developed to predict overall, three dimensional, 
thermal and mechanical properties. This analytical technique was implemented in a user- 
friendly, personal computer-based, windows compatible code called Tex tile Composite 
Analysis for Design (TEXCAD). TEXCAD was used to analyze plain, 5-hamess satin, and 8- 
hamess satin weave composites along with 2-D braided and 2x2, 2-D triaxial braided 
composites. The calculated overall stiffnesses correlated well with available 3-D finite element 
results and test data for both the woven and the braided composites. Parametric studies were 
performed to investigate the effects of yam size on the yam crimp and the overall thermal and 
mechanical constants for plain weave composites. The effects of braid angle were investigated 
for the 2-D braided composites. Finally, the effects of fiber volume fraction on the yam 
undulations and the thermal and mechanical properties of 2x2, 2-D triaxial braided composites 
were also investigated. 

Kev Words; textile composites, modeling, plain weave, satin weave, triaxial braid, yam, 
crimp, design, thermal, stiffness, elastic properties. 


INTRODUCTION 


Composite materials reinforced with woven and braided fabric preforms are being 
considered for potential structural applications in the aircraft and automotive industries. 
Fabric-reinforced 'textile' composites potentially have better out-of-plane stiffness, strength, 
and, toughness properties than tape laminates. They are also amenable to automated and 
efficient manufacturing techniques. However, the architecture of a textile composite is 
complex and, therefore, the parameters controlling its mechanical properties are numerous. 
This necessitates development of analytical models to predict the thermal and mechanical 
properties of textile composites. The present study focuses on the analysis of woven and 
braided textile composites. 

Fabric Architecture 

A two-dimensional (2-D) woven composite laminate is one of the simplest textile 
composites and it consists of stacked, pre-impregnated layers of woven fabric [1] which are 
cured and consolidated by a process similar to tape laminates. The fabric is composed of two 
sets of interlacing, mutually orthogonal (warp and fill) yams (Fig. 1). Each yam is a bundle 
of filaments (or fibers) and the yam size is measured by the number of filaments in the yam. 
The fabric is woven on a loom and its architecture is characterized by the interlacing pattern of 
the warp and fill yams. Some of the commonly used weave architectures are the plain weave, 
5-hamess satin and 8-hamess satin weaves (see Fig. 1). The periodicity of the repeating 
pattern in a woven (or braided) fabric can be used to isolate a small repeating unit cell (RUC) 
which is sufficient to describe the fabric architecture. The RUC for the plain, 5- and 8-hamess 
satin weaves are indicated in Fig. 1 by dark borders. The interlacing of the yams in a fabric 
leads to yam undulation or 'crimp'. Composite stiffness properties are governed by (i) weave 
parameters such as weave architecture, yam sizes, yam spacing, and yam crimp; (ii) laminate 
parameters such as stacking orientations and overall fiber volume fraction; and (iii) material 
parameters such as yam and resin stiffnesses. 
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Braided preforms are fabricated over a cylindrical mandrel. A single layer of a 2-D 
braid consists of either two or three intertwined yarns. The braider yams follow the +0 and 
-6 directions and usually interlace in either alxlora2x2 pattern (see Fig. 2). A 1 x 1 
pattern is similar to the interlacing in a plain weave fabric (see Fig. 2(a)). In a 2 x 2 pattern a 
+0 braider yam continuously passes over two -0 braider yams and then under two -0 braider 
yams and vice versa. A 2-D triaxial braid consists of axial yams in addition to the off-axis 
braider yams (see Fig. 2(b)). The axial yams follow the longitudinal direction and are 
inserted between the braider yams. 

A 2-D braided composite laminate of a certain desired thickness is formed by 
overbraiding layers of the desired architecture on top of each other. After the preform is 
removed from the mandrel, slit along the axial direction, flattened and border stitched, resin is 
introduced by a resin transfer molding process. This leads to a two dimensional material as 
there are no interlayer yams through the laminate thickness. Composite stiffness properties are 
governed by (i) braid parameters such as braid architecture, yam sizes, yam spacing, and yam 
crimp; and (ii) material parameters such as yam and resin stiffnesses and overall fiber volume 
fraction. 

Typical photomicrographs for a 2x2, 2-D braided laminate are shown in Fig. 3. The 
longitudinal section (Fig. 3(a)) along the axial yams indicates minor crimp in these yams. 

Fig. 3(b) shows a transverse section of the laminate. The lenticular cross-sectional shape of 
the axial yams is evident. Also, the axial yams do not stack on top of one another. This is 
due to nesting between layers during the brading operation. Fig. 3(c) shows a section along 
the braider yams. The degree of crimp in the braider yams as they intertwine is well 
illustrated in this section. 

Analytical Techniques 

A number of analytical and numerical techniques have been used to analyze both woven 
and braided composites. All of the proposed models exploit the periodicity of the textile 
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architecture to isolate a representative unit cell (RUC). The RUC is assumed to lie in the 
interior of the composite away from edge effects. Overall properties are calculated by making 
certain assumptions about the RUC internal geometry and deformation field. A review of the 
analytical models used in the past was presented by Raju, Foye, and Awa [2]. Some of the 
earlier models for analyzing woven and braided composites are highlighted below. 

Chou and Ishikawa [3] presented three different models for 2-D woven composites 
based on classical lamination theory (CLT) - the mosaic model, the fiber undulation model, 
and the bridging model. The mosaic model neglected the undulation of the yams and idealized 
the RUC as an assemblage of asymmetric cross-ply laminates. The fiber undulation model 
considered the undulation in the fill yam direction but neglected the undulation in the warp 
yams. The bridging model extended the applicability of the fiber undulation model to the 
analysis of satin weave composites (Fig. 1) by separately modeling the straight portions of the 
yams. 

The fabric geometry model (FGM) [4] was another model based on CLT. It combined 
textile engineering methodology and a modified laminate theory. The geometric model 
consisted of straight yams oriented in six different directions. Such a unit cell representation 
neglected the crimp in the yams. The FGM was also used together with an energy analysis [5] 
to predict stiffness properties of three-dimensional (3-D) braided composites. The FGM model 
was further modified by Pastore and Gowayed [6]. The fiber inclination model (FIM) [7] was 
also based on CLT. It treated the fabric unit cell as an assemblage of inclined unidirectional 
laminae. The FIM also neglected the crimp in the yams. Ramnath [8] presented an approach 
that used a variable yam cross-section geometry model along with admissible displacement and 
stress fileds to determine bounds on the elastic properties of 2-D woven composites. 

Naik and Ganesh [9] and Raju and Wang [10] developed models for woven composites 
that accounted for undulations in both the fill and the warp yams. They computed overall 
stiffness properties by assuming either iso-strain or iso-stress or a combination of these and 
analytically integrating through the volume of the RUC. Due to the tedious analytical 



description of yam cross-sectional shapes and yam paths and the analytical volume 
integrations, this approach is not well suited for the analysis of more complicated architectures 
such as 2-D braids and 2-D triaxial braids. 

Three-dimensional (3-D) finite element models (FEM) [11-15] have also been 
developed for the analysis of plain weave composites. The 3-D FEM are highly computer 
intensive and also require considerable time and effort for model generation. Foye [13] 
developed inhomogeneous finite elements and analyzed sub-cells within the RUC to overcome 
these limitations. He analyzed many textile architectures including plain and satin weaves, 

2-D braids, and 2-D triaxial braids. However, he had to manually calculate the stiffness 
properties, after estimating the material directions and volume fractions, for each element in 
the 3-D model. FEM is, therefore, not well suited for performing parametric studies to 
investigate the sensitivity of overall mechanical properties to architecture parameters. 

Naik, Masters and Ifju [16] presented a simple but accurate geometric modeling and 
analysis technique for 2-D triaxially braided composites. They discretely modeled the yam 
centerline paths within the RUC by assuming sinusoidal undulations and used a yam 
discretization scheme together with a 3-D stress averaging procedure to compute overall 
stiffness properties which correlated well with test data over a wide range of architecture 
parameters. The objective of the present study was to extend this modeling approach to the 
development of a general purpose analysis technique that could (i) discretely model the internal 
yam architecture and (ii) calculate overall thermal and mechanical properties for a variety of 
woven and braided composite materials, such as, plain and satin weaves (Fig. 1), 2-D braids 
and 2-D triaxial braids (Fig. 2) and other textile architectures. Another objective of this effort 
was to develop a technique that was accurate enough to provide good estimates of overall 
properties over a wide range of architecture parameters while being simple enough to be 
performed on a personal computer. 
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ANALYTICAL PROCEDURES 

The analysis of a fabric reinforced composite requires, first, a proper three dimensional 
description of the preform architecture. The geometric modeling of the preform architecture 
was performed, in the present study, by utilizing the periodicity of the textile composite to 
isolate a repeating unit cell (RUC). For each yam within the RUC the yam centerline path 
was then described by assuming a flattened lenticular cross-sectional shape for the interlacing 
yams. This assumption was based on the lenticular yam cross-sectional shapes observed in 
micrographs of the composite (see Fig. 3). Details of the geometric modeling technique and 
its applications to various woven and braided architectures are included in the following 
sections. Having described the fiber architecture, the three dimensional effective stiffnesses 
for the composite were computed by first discretizing each yam in the RUC into yam slices 
and then using the material properties, spatial orientation, and volume fraction of each yam 
slice in a volume averaging technique that assumed an iso-strain state within the RUC. 

Geometric Modeling of Plain Weave Composites 

The RUC for the 2-D plain weave composite is shown in Fig. 4. The sectional view in 
Fig. 4 (section A- A) shows the undulations of a warp yam as it crosses over and under the fill 
yams. Only one layer of the plain weave is shown in the sectional view. The analysis of 
stacked woven and braided layers is discussed in a following section. 

The plain weave composite is usually specified by known quantities such as yam 
spacing, a, in the fill and warp directions, yam filament count, n, for the fill and warp yams, 
yam packing density, p d , filament diameter, d f , and overall fiber volume fraction, V f . These 
known quantities were used to determine the unknown quantities such as yam thicknesses, yam 
cross-sectional areas, yam crimp angle and yam undulating paths which are required to 
discretely model each yam within the RUC. For the sake of simplicity, the yam spacings and 
yam fiber counts for the fill and warp directions were assumed to be equal in the following 
derivation. 
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The yam cross-sectional area, A, was assumed to be the same for both the warp and fill 
yams and also to remain constant along the entire yam path. The projected length, Lp, of the 
yam path for each yam was a function of the yam spacing, a, and was given by Lp = 2 a (see 
Fig. 4). The volume occupied by the four yams within the RUC was given by 4 A Lp. The 
dimensions of the RUC were Lp x Lp x H, where H was the RUC thickness. The volume of 
the RUC not occupied by the yams was assigned to the interstitial matrix. For a yam packing 
density of p d , the overall fiber volume fraction, V f , for the RUC in Fig. 4 can be derived as: 


V, = 


2 P„ A 
Ha 


( 1 ) 


Using the yam filament count, n, and filament diameter, d f , the unknown cross-sectional area, 
A, for the yams was calculated as: 


A 


n d^ n 
4 P„ 


( 2 ) 


For a given overall fiber volume fraction, the unknown thickness, H, was calculated using 
Eqns. (1) and (2). The quantity p d may either be determined by image analysis of 
photomicrographs of sections through the thickness of the composite or may be assumed to be 
between 0.7 - 0.8. The yam thickness, t, was related to the RUC thickness, H, by t = H/2. 
The yam width, w, was calculated by assuming complete coverage (i.e. no gap between 
adjacent yams) and was given by w - a. 

The yam path for each of the four yams in the RUC consists of two straight portions 
and three undulating portions. The undulations are centered around cross-over points (COP) 
where a yam crosses over or under another yam (see Fig. 4). In each of the three undulating 
portions, the yams were assumed to follow a sinusoidal path. The sine function used to 
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describe the yarn centerline path had its origin at the COP and was determined by the vertical 
shift, V g , and the undulating length, L u . For a warp yam, for example, the undulation, Z c , at 
each COP was expressed as: 


Z 


C 



(3) 


where, X c was measured from the corresponding COP along the warp yam direction. A 
negative sign was used in Eqn.(3) to describe the undulation, at the central COP, for the warp 
yam shown in section A-A (see Fig. 4). A positive sign was used in Eqn. (3) to describe the 
undulations at the COPs located on the RUC edges for the same warp yam. Only half the sine 
wave was used to describe the yam centerline paths at the edge COPs. The centerline path for 
the other warp yam in the RUC was described in a similar manner by using Eqn. (3) with the 
appropriate sign and sine wave portion at each COP. The undulations in the fill yams were 
also described by Eqn. (3) by measuring X c along the fill yam direction (i.e. along the Y- 
axis). 

For the plain weave architecture, the vertical shift, V g , in Eqn. (3) was equal to the 
thickness, t, of the yams. The parameter, L u , was determined by considering the cross- 
sectional shape of the yams which consisted of a central flat portion of thickness, t, and two 
sinusoidal lenticular end portions as shown in Fig. 4. The curved portions of the yam cross- 
section were assumed to be described by the sinusoidal functional form of Eqn. (3). Thus, the 
width of the curved portion of the yam cross-section was equal to L u /2. The cross-sectional 
area, A, was expressed as: 


A 


= w t - L u 





(4) 
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The unknown parameter, L u , was calculated using Eqns. (2) and (4). The total length of the 
straight portions, L gt , of each yam was expressed as: 

L st = 2a-2L U (5) 

Thus, Eqns. (1) - (5) were used to determine all the required geometry parameters for the plain 
weave RUC, which included, the yam cross-sectional area, yam thickness, and, yam paths 
based on a knowlegde of the yam filament counts, yam spacing, yam packing density, 
filament diameter, and overall fiber volume fraction. Note that, Eqn. (4) is valid for w > L u 
which physically implies that the parameter, L u , (also the total width of the curved portions of 
the yam cross-section) cannot exceed the spacing, a, between the yams. 

The undulation in the yams is often described by its "crimp angle". In the present 
study, the crimp angle (see Fig. 4) was defined as the positive angle between the tangent to the 
sinusoidal yam centerline path and the X-Y plane (see Fig. 4) at a COP. Thus, the crimp 
angle, 0 C , was expressed as: 


tan(0 c ) 




( 6 ) 


In general, 0 min < 0 C ^ (rc/2), where, 0 min , was the minimum 0 C for a given plain weave 
architecture and was determined by the constraint w > L u . 

Geometric Modeling of Satin Weave Composites 

The RUCs for the 5-hamess satin weave and the 8-hamess satin weave are shown in 
Figs. 5 and 6, respectively. The sectional views in Figs. 5 and 6 (section A-A) show the 
undulations of a warp yam as it crosses over and under the fill yams. 
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The architecture modeling for the satin weave composite was similar to the plain weave 
composite. The yam cross-sectional area, A, was assumed to be the same and to remain 
constant for each yam in the RUC. The projected length, Lp, of the yam path for each yam 
and, A, were used to calculate the volume occupied by the yams within the RUC. The 
volume occupied by the yams within the RUC was given by 2 h A Lp where h = 5 for a 5- 
hamess satin weave and h = 8 for a 8-hamess satin weave. The projected length, Lp, was a 
function of the yam spacing, a, and was given by Lp = h a (see Figs. 5 and 6). The 
dimensions of the RUC were Lp x Lp x H, where H was the RUC thickness. The volume of 
the RUC not occupied by the yams was assigned to the interstitial matrix in the RUC. The 
yam width, w, was related to the yam spacing, a, by assuming complete coverage and was 
given by w = a. For a yam packing density of p d , the overall fiber volume fraction, V f , for 
the 5- and 8-hamess satin weave architectures can be shown to be given by Eqn. (1). The yam 
cross-sectional area, A, for both the satin weave architectures was calculated using Eqn. (2). 
The RUC thickness, H, was then calculated using Eqns. (1) and (2) for both satin weaves. 

Yam thickness, t, was given by t = (H/2). 

The yam undulations were centered around either two or three cross-over points for the 
satin weaves. For example, for the 5-hamess satin weave, the warp yams in the first, third, 
and fourth rows from the top (Fig. 5) had two COPs, while warp yams in the second and fifth 
(section A-A) rows had three COPs. The sine function given by Eqn. (3) was used to describe 
the yam centerline paths at the COPs with the appropriate sign and sine wave portion at each 
COP. 

For both the satin weave architectures, the vertical shift, V g , in Eqn. (3) was equal to 
the thickness, t, of the yams. The parameter, L u , which was also the total width of the curved 
portion of the cross-sectional shape of the yams (see Fig. 4) was calculated using Eqns. (2) 
and (4). The total length of the straight portion of each yam was given by: 

L st = h a - 2 L u (7) 
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where h = 5 for 5-hamess satin weaves and h — 8 for 8-harness satin weaves. Thus, 

Eqns. (1) - (4), (7) were used to determine the geometry of the 5- and 8-hamess satin weave 
composites. The crimp angle, 0 C , was given by Eqn. (6) and the constraint w > L u was used 
to determine 0 min . Note that, Eqn. (7) can be substituted for Eqn. (5) by using h = 2 for the 
plain weave composite, thereby, generalizing the formulation for the plain, 5-hamess satin, 
and, 8-hamess satin weave composites. 

Geometric Modeling of 2-D Braided Composites 

The RUC for the 2-D braided architecture is shown in Fig. 7. The sectional view in 
Fig. 7 (section A-A) shows the undulations of a +0 braider yam as it crosses over and under 
the -0 braider yams. The architecture modeling for the 2-D braided composite was similar to 
the plain weave composite. However, unlike the plain weave, the braider yams, in a 2-D 
braid, interlace at non-orthogonal angles. This fact was accounted for in determining the yam 
paths in a 2-D braid. The yam projected length, Lp, was a function of the yam spacing, a, 
(perpendicular to the yam direction) and the braid angle, 0 br , and was given by: 


L r = 


2 a 


sin(20 br ) 


( 8 ) 


where 0 br was measured with respect to the Y-axis (see Fig. 7). 

As in the plain weave, the yam cross-sectional area, A, was assumed to be the same 
and to remain constant for each yam in the RUC. The projected length, L p , of the yam path 
for each yam and, A, were used to calculate the volume occupied by the yams within the 
RUC. The volume occupied by the four yams within the RUC was given by 4 A Lp. The 
dimensions of the parallelepiped RUC were Lp x Lp x H, where H was the RUC thickness and 
included angle between the two sides with length, Lp, was 20 br . The volume of the RUC not 
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occupied by the yams was assigned to the interstitial matrix in the RUC. The yam width, w, 
(perpendicular to the yam direction) was related to the yam spacing, a, by assuming complete 
coverage and was given by w = a. For a yam packing density of p d , the overall fiber volume 
fraction, V f , for the 2-D braid can be shown to be given by Eqn. (1). The yam cross-sectional 
area. A, for the 2-D braided architecture was calculated using Eqn. (2). The RUC thickness, 
H, was then calculated using Eqns. (1) and (2). Yam thickness, t, was given by t = (H/2). 

As in the plain weave architecture, the yam undulations for the 2-D braid are centered 
around three cross-over points for each yam. The sine function given by 


Z = ± — L sin 
2 


7tX c sin(20 br ) 


V 


(9) 


where, X c is measured from the corresponding COP along the braider yam direction, was used 
to describe the yam centerline paths at the COPs with the appropriate sign and sine wave 
portion at each COP. The vertical shift, V s , in Eqn. (9) was equal to the thickness, t, of the 
braider yams. The cross-sectional shape of the braider yams (perpendicular to the braider 
yam direction) was similar to that shown in Fig. 4 and consisted of a central flat portion of 
thickness, t, and two sinusoidal lenticular end portions of width, L u /2. The total width, L u , 
was calculated using Eqns. (2) and (4). The total length of the straight portion of each yam 
path, L st , (along the braider yam direction) was given by 


= 2 (a - L„) 
sin(20 br ) 


( 10 ) 
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Thus, Eqns. (1), (2), (4), (9) and (10) were used to determine the geometry of the 2-D braided 
composites. The crimp angle, 0 C , was determined by using Eqns. (6) and (9) and the 
limitation w S L u was used to determine 0 min . 

Geometric Modeling of 2-D Triaxial Braided Composites 

The RUC used for the 2x2, 2-D triaxial braid is shown in Fig. 8. The sectional view 
A-A shows the undulation of a braider yam which undulates over and under the axial yams in 
the RUC. The 2-D triaxial braided composite is usually specified by known quantities such as 
braid angle, 0 br , axial yam spacing, d 0 , yam filament counts for the axial and braider yams, 
n Q and n b , respectively, yam packing density, p d , filament diameter, d f , and overall fiber 
volume fraction, V f . These known quantities were used to determine the unknown quantities 
such as yam thicknesses, yam cross-sectional areas, yam crimp angle and yam undulating 
paths which are required to discretely model each yam within the RUC. A detailed description 
of the geometric modeling approach for the 2x2, 2-D triaxial braided composite was presented 
in Ref. 16. Simplified equations for the same geometric model are presented here. 

As before, the yam cross-sectional areas, A 0 , and A b , for the axial and braider yams, 
respectively, were assumed to remain constant along the entire yam path. The projected 
lengths, L 0 and 1^, of the axial and braider yam paths were a function of the axial yam 
spacing, d 0 , and the braid angle, 0 br , and were given by L c = 8 d 0 cot(0 br ) and 
Lj, = 2 d o /sin(0 br ), respectively. Note that L 0 represents the total projected length of the four 
axial yams within the RUC. The total volume occupied by the axial and braider yams was 
given by (A 0 L 0 + 8 Lj, A b ). The dimensions of the parallelepiped RUC were L*, x 1^, x H, 
where H was the RUC thickness and the included angle between the two sides with length, Lj,, 
was 20 br . The volume of the RUC not occupied by the yams was assigned to the interstitial 
matrix in the RUC. For a yam packing density of p d , the overall fiber volume fraction, V f , 
for the RUC in Fig. 8 was derived as: 
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( 11 ) 


v = p„ (A 0 sin (28^ ) + 4 A b sinCe^)) 
f Hd o sin(20 br ) 

where, H is the thickness of the RUC. Using the known yam filament counts, n 0 and n b , 
filament diameter, d f , the unknown cross-sectional areas, A c and A b , for the axial and braider 
yams, respectively, were calculated using Eqn. (2). For a given overall fiber volume fraction, 
the unknown thickness, H, was then calculated using Eqn. (1 1). The braider yam thickness, 
th, was a function of the RUC thickness, H, and the yam thickness ratio, m 


t b 


H 

2 + m’ 



( 12 ) 


where ^ was the thickness of the axial yams. The unknown yam thickness ratio, m, was 
determined recursively as described below. The braider yam width, w b , (perpendicular to the 
yam direction) was calculated by assuming complete coverage (i.e. no gap between adjacent 
yams) and was given by 


w b = d 0 cos(0 br ) 


(13) 


The yam path for the axial yams was assumed to be straight throughout the RUC. The 
braider yam path consisted of straight portions and undulating portions. In the undulating 
portion the braider yam was assumed to follow a sinusoidal path at each COP. The sine 
function used to describe the braider yam centerline path, Z c , at each COP was expressed as: 


±±J nX ' s^w)) 

2 l L- J 


(14) 
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where, X c was measured from the corresponding COP along the braider yam direction. 

Eqn. (14) was used to describe the yam centerline paths at the COPs with the appropriate sign 
and sine wave portion at each COP. The vertical shift, V g , in Eqn. (14) was equal to the sum 
of the axial and braider yam thicknesses and was given by V g = ^ + V 

The parameter, L u , was determined by considering the cross-sectional areas, A 0 and 
A b , of the axial and braider yams which were expressed by Eqns. (15) and (16), respectively, 
as: 


A„ = w e t 0 - L u V. (l - ~(s + cos(s))j. 



(15) 


where, w Q = d 0 - 2 s (L u /rt) was the width of the axial yams, and 

A„ = w b t b - - L u (t b s - V. (1 - cos(s))) (16) 

71 

The unknown parameters, L u and m, were determined by a recursive solution of Eqns. (15) 
and (16) in which A 0 and A b were calculated using Eqn. (2). The total length of the straight 
portion of each braider yam was given by 


2 (d, - L u ) 
sin(0j 


(17) 


Thus, Eqns. (2), (11) - (17) were used to determine the geometry of the 2x2, 2-D, triaxial 
braided composites. The crimp angle, 0 C , was given by Eqns. (6) and (14) and the limitations 
w b > L u (2 s/re) and w 0 > L u (1 - (2 s/re)) were used to determine the 0 min . 
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Discretization of Yarns 

Overall composite properties were determined by discretizing each of the yams within 
the RUC. The straight portions of each yam path were modeled as a single slice with length, 
L rt . Along an undulating portion, the yam was divided into n, equal, piecewise straight slices 
made perpendicular to its in-plane direction and normal to the X-Y plane (Fig. 8). Thus, the 
undulating, sinusoidal yam path was approximated by n interconnected straight yam slices. 
The volume of each yam slice was computed as (A y Ly/n) where A y and L y are the yam cross- 
sectional area and the length of the undulating portion, respectively. Note that, L y = L u for 
the woven architectures (Figs. 4-6) and L y = L u for the braided architectures (Figs. 7, 8). 

The spatial orientation of each yam slice was described by the in-plane angle, 0, that it 
made with the X-axis (Fig. 8) and the out-of-plane angle, P, that it made with the X-Y plane. 
For the woven architectures, the angle 0 was either 0 (warp yams) or 90 (fill yams) degrees. 
For the braided composites, the angle 0, for the braider yams, was determined from a 
knowledge of the braid angle. The angle P, for both woven and braided architectures, was 
calculated, for each yam slice, by differentiating the sine function used to describe the 
undulating yam centerline path. For the yam slice that modeled the straight portions of the 
yam path, P = 0 was used. 

Thus, each yam in the RUC was approximated by straight yam slices with known 
volumes and known orientation angles 0 and p. The total volume occupied by all the yam 
slices was subtracted from the volume of the RUC to compute the volume occupied by the 
interstitial matrix pockets in the RUC. The interstitial matrix was represented as an isotropic 
material slice having orientation angles 0 = P = 0. 

Calculation of Three-Dimensional Effective Stiffnesses 

For the purpose of calculating overall properties, the RUC of the textile composite was 
treated as a spatially oriented fiber composite composed of yam slices with transversely 
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isotropic material properties and longitudinal material axes oriented at known angles 0 and 0. 
A stress averaging technique [17] based on an iso-strain assumption within the RUC was used 
to compute overall composite properties as described below. 

The three-dimensional thermoelastic stress-strain relations, in global XYZ coordinates, 
for the composite RUC were defined using effective material properties to relate the average 
values of the state variables. Using the usual contracted notation [18], the stress-strain relation 
was written in matrix form as 

{a} = [C eff ] ({i} - {a} AT) (18) 

where {a} = {a„ , , o n , , a lx }, the overbar denotes average values over the 

RUC, {i} = {i^, i^, i^, e^, i^, i^}, {a} = {a^, a^, a^, 0, 0, o} are the 

overall coefficients of thermal expansion (CTE), AT is the change in temperature and [C g ^] 

is a 6x6 matrix [18] of the overall, three-dimensional effective stiffnesses of the composite. 

For a set of displacement boundary conditions applied to the RUC the average strains 
in the composite were given by {e}. The following procedure was used to determine {c}. 

An iso-strain state was assumed to exist within the RUC. This provided for compatibility 
among yam slices and also led to displacements that satisfied the boundary conditions. In the 
global RUC coordinates (XYZ coordinates), the strains in each yam slice were, therefore, 
given by 



( 19 ) 
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where the subscript m denotes the yam slice. Knowing the strains and assuming an iso- 
thermal state in the RUC, yam slice stresses in the local yam slice material coordinate system 
(123 coordinates) (see Fig. 8) were computed as 

{O }. = [c ]. ({e }. - {«'}. at) 

where, (20) 

{e L = [T] m {e}„ 

The 6x6 stiffness matrix [c’j defines the three-dimensional stress-strain relationship for the 
m* yam slice which is assumed to be a transversely isotropic material and (a } are the CTE 

m 

for the m** 1 yam slice. The five independent material constants (E n , E^, G 12 , v 12 , and v^, 
subscript 1 corresponds to the axial fiber direction) that are required to define the [c ] m matrix 

and the yam slice CTE a n , ot 22 and a 33 could be measured quantities based on yam tests. 
They could also be estimated by a rule of mixtures or a micromechanics analysis from 
constituent properties. Such an estimate would require a knowledge of fiber properties, matrix 
properties and the yam packing density p d (or yam fiber volume fraction). The transformation 
matrix [T] m for the strains in Eq. (20) is defined in Appendix A. 

The average stresses in the RUC were found by first transforming yam slice stresses to 
the global RUC coordinates and then computing the volume average of the stresses in all the 
yam slices as 

{a} = £v m {o} m 

m=l 

where, (21) 

M„ = ml {o'}. 
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where V m is the volume fraction of the m 01 yam slice in the RUC and N is the total number of 
yam slices (including the interstitial matrix material slice) in the RUC. The [T]^ matrix is the 

transpose of the matrix [T] m defined in Eqn. (Al) and is used to transform stresses from yam 
slice material coordinates (123) to the global RUC coordinates (XYZ). 

Based on Eqns. (18) - (21) the effective stiffness matrix [C^.] of the RUC was written 
in terms of the yam slice stiffness matrices, transformation matrices and volume fractions as 

[C eff ] - £(vjT]l[c']jT].) (22) 

and the overall CTEs for the RUC were given by 



{«'} )} 

^ m=l\ L Jm m / J 


(23) 


A convergence study was performed to determine the appropriate number of yam 
slices, n, in the undulating portions of the yams, required to yield converged values for the 

overall stiffnesses [C g ^]. The computed stiffnesses for the woven and braided composites in 
the present study were found to be unchanged for n £ 12. A value of n = 12 was used for all 
the analyses in this study. The overall stiffness matrix [C e ff ] was inverted to obtain the 

overall compliance matrix [S g ^ ] which was used to determine overall moduli and Poisson's 
ratios [18]. 

Analysis of Stacked Oriented Textile Layers 

For laminates made by stacking layers with different orientations, the thermo-elastic 
properties for each layer can be obtained from the models described in the preceeding sections, 
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using appropriate tensor transformations. For the layered textile composite laminate, it is 
assumed that classical lamination theory can be used to compute thermal and mechanical 
laminate properties. Each layer of the textile is treated as a lamina in this analysis. Nesting 
between the layers, by which layers slide into the available space in the adjacent layers, is not 
accounted for in the present framework of the classical lamination scheme. 

RESULTS AND DISCUSSION 

The above geometric modeling and analytical technique, for textile composite overall 
thermo-mechanical properties, was implemented in a user-friendly, windows compatible, 
personal computer-based code called Textile Composite Analysis for Design (TEXCAD). 
TEXCAD is also capable of predicting textile composite damage onset, damage growth, failure 
modes and ultimate strength. The formulation for the strength analysis will be presented in a 
subsequent paper. In the present paper, TEXCAD was used to analyze 2-D woven (plain, 5- 
and 8-harness satin), 2-D braided, and 2x2, 2-D triaxial braided composites. The calculated 
overall stiffnesses were correlated with available 3-D finite element (FEM) results and test 
data. Parametric studies were also performed to study the effects of various fabric parameters 
on overall properties. 

All the analyses in the present study were made for yams made of Hercules AS4 
graphite fibers impregnated with Hercules 3501-6 epoxy matrix. A yam fiber packing density, 
p d , of 0.75 and a fiber diameter, d f , of 0.007 mm was used in all calculations. The thermo- 
mechanical properties used for the impregnated yams (V f = 75%) and the resin are given in 
Table 1. The CTE's for the impregnated yams with a 75% fiber content were calculated using 
a micromechanics analysis based on an Airy's stress function approach for a square fiber 
packing [21]. The CTE for the 3501-6 epoxy resin (40x1 O^C) was available in the literature 
[22]. The CTE's for the fiber (<x u = -0.61xl0 6 /°C and 042 = 3.8xl0^/°C) were back- 
calculated using the square array micromechanics model [21] together with the available 
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CTE’s for the resin and for a unidirectional 60% fiber volume fraction AS4/3501-6 laminate. 
The laminate CTE's were taken from manufacturers data sheets [23]. 

Woven Composites 

The thermal and mechanical elastic constants for plain, 5-hamess satin and 8-harness 
satin weave composites calculated using TEXCAD (the present analysis) are shown in Table 2. 
The calculated results correlated well with finite element (FEM) results and test data reported 
in Ref. 13. A yam spacing of 1.41 1 mm and an overall fiber volume fraction, V f , of 64% 

[13] were used for all the three woven architectures. The yam size was not specified in Ref. 
13. It was, therefore, assumed to be 3 k (k - one thousand filaments) for both the fill and 
warp yams. Using the geometric modeling approach described earlier, a 3 k yam size led to a 
19.4 degree crimp angle for the given yam spacing and V f . The crimp angle, for the plain 
weave composite, measured from the micrographs in Ref. 13, varied between 8 and 18 
degrees. A 3 k yam size was also used for calculating the results in Table 2 for the 5-hamess 
satin and 8-hamess satin weave composites. 

Table 2 also includes results for a [0/90/r] s laminate (r - resin layer) calculated using 
classical lamination theory (CLT). To allow a direct comparison with the woven composites, 
the 0 and 90 degree plies were assigned the impregnated yam (V f = 75%) material properties 
(Table 1). Since the yams (plies) had a 75% fiber content, a resin layer, r, was added to the 
cross-ply laminate to achieve an overall composite V f = 64%. The volume fraction of the 
resin layer in the laminate was equal to that of the interstitial matrix in the woven composites. 
The CLT results further verify the trends in the results for the three architectures. Since the 
plain, 5-hamess satin and 8-hamess satin weaves have increasing amounts of yam straight 
sections, L st , (see Eqns. (5) and (7)) their moduli and CTE's should approach those of the 
cross-ply laminate. 

Woven composite prepreg is often manufactured in different yam sizes, such as, 3 k, 

6 k, and 12 k. A parametric study was performed, using TEXCAD, to investigate the effects 
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of yam size on the overall thermal and mechanical properties. Figure 9 shows the effects of 
yam size on the crimp angle and the stiffness properties for a plain weave composite. For the 
same yam spacing, yam packing density and overall composite volume fraction, the crimp 
angle increased significantly with increasing yam size. The in-plane moduli, E^ and Ey y , 
decreased while, E^, increased with increasing yam size. The shear moduli, the in-plane 
Poisson's ratio, v xy , (Fig. 9) and the in-plane CTE's, and a^, (see Fig. 10) were not very 
sensitive to yam size effects. However, the through-thickness CTE, a^, (Fig. 10) and 
Poisson's ratios, v K and v yz , (Fig. 9) were very sensitive to yam size. 

Note that, since the overall composite V f and interstitial matrix volume fraction was 
kept constant while the yam size was varied, a CLT-based analysis would have predicted no 
change in the overall properties with varying yam size. It is, therefore, necessary to explicitly 
model the yam architecture and crimp to investigate the effects of fabric parameters on overall 
thermal and mechanical properties of textile composites. 

2-D Braided Composites 

A yam spacing of 1.41 1 mm (perpendicular to the yam axial direction) and an overall 
fiber volume fraction, V f , of 64% were used for the 2-D braided composites analyzed. As in 
the plain weave composite, a 2 k yam size, a 0.75 yam packing density, and a 0.007 mm fiber 
diameter were used for calculating the results for the different 2-D braided composites. A 
parametric study was performed to study the effects of braid angle variations on composite 
properties. Figures 1 1 and 12 show the effects of braid angle on the crimp angle, thermal, and 
mechanical properties of the 2-D braided composites. A good correlation was obtained with 
FEM results [13] for a 45 degree braid angle for the same braid parameters (Fig. 11). Since 
the 45 degree braided composite is equivalent to a plain weave composite rotated by 45 
degrees (see Figs. 4 and 7), the computed results for a 45 degree braided composite were also 
verified by transforming the plain weave stiffnesses by 45 degrees. 
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The in-plane properties, E^, Ey y , G xy and v xy , varied significantly with braid angle 
(Fig. 11). Through-thickness Poisson's ratios, v M and v yz , and modulus, E^, were less 
sensitive to braid angle. Crimp angle (Fig. 12) varied between 6.7 and 13.2 degrees for braid 
angles between 15 and 75 degrees. The in-plane CTE’s, a**, and a yy , varied significantly 
with braid angle. Through-thickness CTE, a u , was less sensitive to braid angle. 

2x2, 2-D Triaxial Braided Composites 

The three different triaxial braided composites tested in Ref. 19 were analyzed using 
the present TEXCAD analysis. The yarn architecture parameters such as braid angle, braider 
and axial yarn sizes, axial yarn spacing and overall composite volume fraction for the three 
architectures (Al, Bl, B2) are given in Table 3. These parameters along with a yarn packing 
density of 0.75 and a fiber diameter of 0.007 mm were used as input to TEXCAD for 
calculating the results in Table 4. The specimens in Ref. 19 were made with AS4 graphite 
fiber yams and Shell 1895 epoxy resin. It was assumed that the 1895 epoxy had the same 
material properties as the 3501-6 resin (Table 1). Thus, the material properties for the 
impregnated yams and the resin in Table 1 were used in the present analysis of 2-D triaxial 
braids. 

Correlations of the computed elastic constants (Table 4) with test data [19] and FEM 
results [13, 19] were very good for all the three architectures. Further testing will be required 
to verify the results for the in-plane shear moduli, the through-thickness properties and the 
CTE's. It was shown in Ref. 16 that the results for in-plane moduli (including G xy ) and 
Poisson' ratio using the present analysis technique correlated well with test data over a wide 
range of architecture parameters for 2x2, 2-D triaxial braided composites. Parametric studies 
were also conducted in Ref. 16 to investigate the effects of braid angle, yam size, and axial 
yam content on overall stiffness properties of triaxial braided composites. 

The overall fiber volume fraction for a textile composite, measured after fabrication, 
can often differ by a few percent from the specified fiber content. Furthermore, variations 
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may occur in the fiber volume fraction from one composite plate to another. It is, therefore, 
important to investigate the sensitivity of overall thermal and mechanical properties to fiber 
volume fraction. Figures 13 and 14 show results (for the B1 architecture) for the effects of 
fiber volume fraction on the thermal and mechanical properties of triaxial braided composites. 
Except for and Ey y , which increased proportionately with V f , the shear moduli, E^, the 
Poisson's ratios and the in-plane CTE's (Fig. 14) were not very sensitive to variations in V f . 
For a 9% increase in V f (55% to 60%), and Ey y increased by about 8.3%. This trend in 
the Ejm was further verified by test data (see Fig. 13) for the B1 architecture [20], Figure 14 
shows that the crimp angle increased while the a n decreased significantly with increasing V f . 

SUMMARY 

A general purpose micromechanics analysis that discretely models the yam architecture 
within the textile repeating unit cell, was developed to predict overall, three dimensional, 
thermal and mechanical properties. This analytical technique was implemented in a user- 
friendly, personal computer-based, windows compatible code called Tex tile Composite 
Analysis for Design (TEXCAD). TEXCAD was used to analyze plain, 5-hamess satin, and 8- 
hamess satin weave composites along with 2-D braided and 2x2, 2-D triaxial braided 
composites. The calculated overall stiffnesses correlated well with available 3-D finite element 
results and test data. 

Based on a parametric study performed on plain weave composites, it was found that 
for the same yam spacing, yam packing density and overall composite volume fraction, the 
crimp angle increased significantly with increasing yam size. The in-plane moduli, E xx and 
Ey y , decreased while, E^, increased with increasing yam size. The shear moduli, the in-plane 
Poisson’s ratio, v xy , and the in-plane CTE's, a xx and a yy , were not very sensitive to yam size 
effects. However, the through-thickness CTE, a zz , and Poisson's ratios, v K and v yz , were 
very sensitive to yam size. 
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The effects of braid angle were investigated for the 2-D braided composites. The in- 
plane properties, E^, Ey y , G xy and v xy , varied significantly with braid angle. Through- 
thickness Poisson's ratios, v u and v yz , and modulus, were less sensitive to braid angle. 
Crimp angle varied between 6.7 and 13.2 degrees for braid angles between 15 and 75 degrees. 
The in-plane CTE's, a H , and Oy y , varied significantly with braid angle. Through-thickness 
CTE, a^, was less sensitive to braid angle. 

The effects of fiber volume fraction on the thermal and mechanical properties of 2x2, 
2-D triaxial braided composites were also investigated. Except for E^ x and Ey y , which 
increase proportionately with V f , the shear moduli, E^, the Poisson's ratios and the in-plane 
CTE’s were not very sensitive to variations in V f . For a 9% increase in V f (55% to 60%), E^ 
and Ey y increased by about 8.3%. The crimp angle increased while the decreased 
significantly with increasing V f . 

The present TEXCAD analysis provides a useful tool for the prediction of thermo- 
elastic constants for a wide range of textile composite architectures. It also provides a means 
to conduct sensitivity studies for the numerous parameters that can affect the properties of 
textile composites. Finally, it has the potential for being implemented as an element 
constitutive model in a finite element formulation, thus, providing the capability of structural 
analysis for textile composite structures. 
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APPENDIX A 


The transformation matrix [T] m for the strains in Eq. (20) is defined as [16] 


*n 


a I3 

a i2 a 13 

a ii a t3 

a n a 12 

*2. 

a 2 * 

*L 

a 22 a 23 

a 23 a 21 

321*22 

*31 

a 32 

a 33 

a 32 a 33 

a 33 a 31 

a 31 a 32 

2a 2 ,a 31 

2a 32 a 22 


( a 22 a 33 +a 23 a 3 2 ) 

( a 23 a 3I a 2i a 33 ) 

(*21*32 ■*" *22*31 ) 

2a u a 3 , 

2 a i2 a 32 

2 a i3 a 33 

( a 32 3 13 a 33 a !2 ) 

( a ii a 33 a !3 a 31 ) 

(*31*12 *32*11 ) 

2a,,a 2I 

^ a i2 a 22 

23)3*23 

( a !2 a 23 ^ a i3 a 22 ) 

( a !3 a 2I + a ii a 23 ) 

(*11*22 ^ *12*21 ) 


(Al) 

where a^ (i = 1-3, j = 1-3) are the direction cosines between the yam slice material 
coordinate axes (123) and the global RUC axes (XYZ) (see Fig. 8). For a yam slice with its 
longitudinal material axis (1-axis) oriented at angles 0 and P with respect to the global RUC 
coordinates, the 3x3, [a^] matrix (i = 1-3, j = 1-3) is given by 



cos(0)cos(P) sin(0)cos(P) sin(P) 
-sin(0) cos(0) 0 

-cos(0)sin(P) -sin(0)sin(p) cos(P) 


(A2) 


Equation (A2) was derived by assuming that the 2-axis (i.e. the in-plane transverse direction, 
see Fig. 8) of the yam slice was always perpendicular to the global Z-axis for any 0 or P 
rotation of the yam slice 1-axis (i.e. the axial fiber direction). This merely simplified the 
definition of the direction cosines without compromising the generality of the approach. 
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Table 1.- Yam and resin properties used in the analysis [13, 16, 21, 22, 23]. 


Material 

E n , GPa 

E 22 > GPa 

Gj 2 > GPa 

v 12 

v 23 

CTE a u , 
10</°C 

CTE 022, 
lO^C 

Yam 

144.80 

11.73 

5.52 

0.23 

0.30 

-0.324 

14.000 

Resin 

3.45 

3.45 

1.28 

0.35 

0.35 

40.000 

40.000 


Note: The 1-axis is along the fiber direction, the 2-axis is perpendicular to the fibers but in 
the plane of the lamina and the 3-axis is in the out-of-plane direction. 


Table 2.- Comparison of results for woven composites. 


Laminate 

type 

Approach 


E zz> 

GPA 

G xz> 

G yz> 

GPa 

G xy> 

GPa 

v 

v xz> 

^yz 

v *y 

®xx» ®yy> 

lO-W 

10^/°C 

Plain 

TEXCAD 

64.38 

11.49 

5.64 

Bl 

.396 


1.334 

20.712 

weave 

FEM [13] 

63.78 

11.38 

4.97 

m 

.329 

Mm 

- 

- 


Test [13] 

61.92 

- 

- 

■ 

- 

.110 

- 

- 

5-hamess 

TEXCAD 

66.33 

11.51 

4.93 

4.89 

.342 


1.462 

21.242 

satin 

FEM [13] 

65.99 

11.38 

5.03 

4.96 

.320 

1 

. 


weave 

69.43 

- 

- 

5.24 

- 

.060 

_ 

_ 

Test [13] 

8-hamess 

TEXCAD 

66.81 

11.51 

1 

4.89 

.329 


1.494 

21.381 

satin 

FEM [13] 

66.74 



4.96 

.320 

1 



weave 

72.19 


Hi 

6.76 

- 

.060 

_ 

_ 

Test [13] 

[0/90/ri], 

laminate 

CLT 

67.60 

- 

- 

4.89 

- 

.038 

1.548 

- 


•The r represents a resin layer used in the laminate to simulate the effect of the interstitial 
matrix in the woven composite. 
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Table 3.- Input parameters for 2x2, 2-D triaxial braided composites [19]. 


2-D triaxial 
braid 

Braid angle, 
degrees 

Braider yam 
size, k 

Axial yam 
size, k 

Axial yam 
spacing, mm 

Overall fiber 
volume 
fraction, % 

A1 

62.3 

12 

24 

6.096 

54.0 

B1 

67.4 

6 

18 

5.320 

51.2 

B2 

67.9 

6 

18 

5.820 

52.0 


Table 4.- Comparison of results for 2x2, 2-D triaxial braided composites. 


Braid 

type 

Approach 

^xx> 

GPa 

Eyy- 

GPA 

Ezz, 

GPa 

G X y, 

GPa 

G xz< 

GPa 

G yz> 

GPa 

v *y 

<*xx> 

lO-W 

CLyy, 

lO^C 

A1 

TEXCAD 

45.14 


10.42 

14.84 

4.64 

4.13 

.291 

1.675 

2.650 

FEM [19] 

53.16 

46.95 

- 

- 

- 


.305 

- 

- 


Test [19] 

45.44 

45.64 

- 

- 

- 

- 

.307 

- 

- 

B1 

TEXCAD 

48.42 

42.48 

10.13 

11.07 

4.40 

3.90 

.205 

1.937 

2.901 

FEM [19] 

62.12 

47.51 

- 

- 

- 

- 

.211 

- 

- 


Test [19] 

44.47 

45.16 

- 

- 

- 

- 

.199 

. 


B2 

TEXCAD 

50.42 

42.63 

10.23 

11.08 

4.31 

3.90 

pS3| 

1.870 

2.897 

FEM 1 [19] 

51.09 

43.58 

12.07 

12.41 

- 

- 

1 

- 

- 


Test [19] 

48.47 

43.51 

- 

- 

- 

- 

.190 

- 

- 


1 The E^ and G xy were reported in Ref. [13] for the B2 architecture. 
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Fig. 1.- Typical yam interlacing patterns for woven fabrics. 
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Fig. 2.- Typical yam interlacing patterns for 2-D braided fabrics. 
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Fig. 5.- The RUC geometry and notation for a 5-hamess satin weave composite. 
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- Cross-over points 
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Fig. 7.- The RUC geometry and notation for a 2-D braided composite. 
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Fig. 9 - Variation of plain weave crimp angle, moduli and Poisson's ratios with 



100 


seaiBep 'd|6ue dujjjQ 



41 


Yam size, k 

Fig. 10.- Variation of plain weave crimp angle and CTE*s with yam size. 
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Braid angle, degrees 

Fig. 12.- Variation of 2-D braided composite crimp angle and CTE's with braid angle. 
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Fig. 14.- Variation of crimp angle and CTEs with V f for 2-D triaxial braided composite. 
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